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Superconductivity in the two component model of coexisting local elec- 
tron pairs (hard-core charged bosons) and itinerant fermions coupled via 
charge exchange mechanism is discussed. The cases of isotropic s-wave and 
anisotropic pairing of extended s-wave and d x 2_ y 2 symmetries are analyzed 
for a 2D square lattice within the BCS-mean field approximation and the 
Kosterlitz-Thouless theory. The phase diagrams and superconducting char- 
acteristics of this induced pairing model as a function of the position of the 
local pair (LP) level and the total carrier concentration are determined. The 
model exhibits several types of interesting crossovers between the BCS like 
behavior and that of LP's. In addition, the Uemura plots are obtained for 
extended s and d x 2_ y 2 pairing symmetries. Finally, we analyze the pairing 
fluctuation effects (in 3D) within a generalized T-matrix approach. Some of 
our results are discussed in connection with a two-component scenario of pre- 
formed pairs and unpaired electrons for high temperature superconductors. 



I. INTRODUCTION 

A mixture of interacting charged bosons (local electron pairs with q = 2e) and electrons 
can show features which are intermediate between those of local pair (bipolaronic) supercon- 
ductors and those of classical BCS systems. Such a two component (boson-fermion) model is 
of relevance for high temperature superconductors (HTS) and other exotic superconductors 
[1-13]. A similar model has also been adopted for the description of a resonance s-wave 
superfluidity in Fermi atomic gases with a Feshbach resonance [14,15]. 



1 



Recently, we have studied a generalization of this model to the case of anisotropic pairing 
[11-13]. Here, we briefly outline the study and present some further results concerning the 
phase diagrams and superconducting properties of such a system in the case of isotropic 
s-wave and anisotropic rf-wave pairing, for a 2D square lattice (Sec.II-III), as well as for 
s-wave pairing for a 3D simple cubic (sc) lattice (Sec. IV). 

II. THE MODEL 

We consider the model of coexisting electron pairs (hard-core bosons "b") and itinerant 
" c" electrons defined by the following effective Hamiltonian 

h = j> k - n)<L<** + 2 E( A o - rit>fo - E J iAh + 

kcr i ij 

E [^q(k) C k+q/2,T C -k+q/2a & q + k - C ] + H d ( X ) 
k,q 

where ek refers to the band energy of the c-electrons, Ao measures the relative position 
of the LP level with respect to the bottom of the c-electron band, \i is the chemical po- 
tential which ensures that the total number of particles in the system is constant, i.e. 
n = jj (X)ko-( c L- c ko-) +2Z)z(&!&i)) = n c + 2ub- n c is the concentration of c-electrons, % 
is the number of local pairs per site. Jy is the pair hopping integral. He denotes Coulomb 
interaction terms. The operators for local pairs b\,bi obey the Pauli spin 1/2 commuta- 
tion rules. Vq(k) describes the coupling between the two subsystems. We will consider the 
case Vq(k) = Vo(k) = Ifa/y/N, and neglect its q dependence at small q. The interaction 
term takes the form of coupling, via the center of mass momenta q, of the singlet pair of 
c-electrons and the hard-core boson 6 q : 

H 1 = ^=Y.^ B \K + b\ B ^)- (2) 
v jy q 

-Sq = Sk 0fe c k+ q /2 t c -k+q/2 j. denotes the singlet pair creation operator of c-electrons and I 
is the coupling constant. The pairing symmetry, on a 2D square lattice, is determined by 
the form of fa, which is constant (1) for on-site pairing (s), <pk = Ik = cos{k x ) + cos(k y ) 
for extended s-wave (s*) and <pk = Vk = cos(k x ) — cos(k y ) for d x 2_ y 2-wave pairing (d). In 
general, one can consider a decomposition I fa = go + g s 1k + 9dVk, with appropriate coupling 
parameters for different symmetry channels. 
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The superconducting state of the model is characterized by two order parameters: xq = 
Y,k 0fc( c lt c -fej.) an d Po = w + 111 the BCS-mean-field approximation (MFA) the 
free energy of the system (for H c = 0) is evaluated to be: 



F,N = ~JN ^ ln [2 cosh (^/2)] - ^ In [2 cosh^A)] + C, 
C=-e b + A + /i(n c + 2n B ) - 2/i - 2I\x \p x + J (p x ) 2 , 



(3) 
(4) 



where the quasiparticle energy of the c-electron subsystem is given by 



£ k = ^el + Al and e k = e k - p, A 2 = / 2 0i(pg) 2 , A = ^(A - p) 2 + (-/|rr | + J pg) 2 - 

— J2i^jJij- ft — l/k B T. For the 2D square lattice the c-electron dispersion is 
e k = e k — e b = — 2i [cos^) + cos(k y )] — 4t 2 cos(k x ) cos(k y ) — €(,, with the nn and nnn 
hopping parameters t and t 2 , respectively, e& = mme k . It should be noted that the energy 
gap in the c-band is due to nonzero Bose condensate amplitude (|(6)| ^ 0), and well defined 
Bogoliubov quasiparticles can exist in the superconducting phase. The order parameters 
and the chemical potential are given by 



dF 



= 0, 



dF 



OF 

°' 7r = - 

op 



(5) 



dx dp x 

The superfluid stiffness derived within the linear response method and BCS theory, for the 
case Jij = 0, is of the form: 
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Ps = 



2N 2 ^\\dk x ) <)F k 2dkl L 



e k . , (3E k 

—— tanh 

2 , 



(6) 



where f(E k ) = 1/ [exp((3E k ) + 1] is the Fermi-Dirac distribution function. In the local limit: 
A~ 2 oc (167re 2 /h 2 c 2 )p s , where A is the London penetration depth. 

The mean-field transition temperature (T c Mi?A ), at which the gap amplitude vanishes, yields 
an estimation of the c-electron pair formation temperature [12,13] and is given by 



1 = 



P tanh(/3, 



N 



c e k 



2cl 



tanh 



™(A - p) 



2(A -/i) 



(7) 



Due to the fluctuation effects the superconducting phase transition will occur at a critical 
temperature lower than that given by the BCS-MFA theory. In 2D, T c can be derived within 
the Kosterlitz-Thouless (KT) theory for 2D superfluids [16], which describes the transition 
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in terms of vortex-antivortex pair unbinding. We evaluate T c using the KT relation for the 
universal jump of the (in-plane) superfluid density p s at T c [16] : 

-k B T c = Ps (T c ), (8) 

71 

where p s (T) is given by Eq.(6) and x (T), pg(T), /i(T) are given by Eqs.(5). Thus, the 
critical temperature denoted further by T^ T is determined from the set of four self-consistent 
equations. In the weak coupling limit (\I \/2D < 1, J = 0, D-the half-bandwidth, / = 
-|/ |), T C KT /T C MFA -> 1 if \I \/2D -> 0. 



III. RESULTS FOR 2D ELECTRON SPECTRUM 

A comprehensive analysis of the phase diagrams and superfluid properties of the model 
Eq.(l) for different pairing symmetries including s, the extended s (s*) and d x 2_ y 2-wave 
symmetries was performed in Refs. [11-13]. Below we will discuss these results including 
some additional ones. (The term He will not be considered). 

In the absence of interactions, depending on the relative concentration of " c" electrons 
and LP's we distinguish three essentially different physical situations. For n < 2 it will be: 

(i) A < such that at T = OK all the available electrons form local pairs (2n B 3> n c ) (LP); 

(ii) A > such that the "c" electron band is filled up to the Fermi level \x = A and the 
remaining electrons are in the form of local pairs (the "c+b" or Mixed regime, < 2n B , n c < 
2) (LP+E); 

(iii) A > such that the Fermi level \x < A and consequently at T = OK all the available 
electrons occupy the "c" electron states (the c-regime or"BCS", n c 3> 2n B ) (E). 

For |J | 7^ 0, in the case (ii) superconductivity is due to the interchange between local 
pairs and pairs of "c" electrons. In this process "c" electrons become "polarized" into 
Cooper pairs and local pairs increase their mobility by decaying into "c" electron pairs. 
In this intermediate case neither the standard BCS picture nor the picture of local pairs 
applies and superconductivity has a "mixed" character. The system shows features which 
are intermediate between the BCS and preformed local pair regime. This concerns the 
energy gap in the single-electron excitation spectrum (£^(0)), the k B T c / E g (0) ratio, the 
critical fields, the Ginzburg ratio k, the width of the critical regime as well as the normal 
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state properties. In case (i) the local pairs can move via a mechanism of virtual excitations 
into empty c-electrons states. Such a mechanism gives rise to the long range hopping of 
LP's (in analogy to the RKKY interaction for s-d mechanism in the magnetic equivalent). 
The superconducting properties are analogous to those of a pure local pair (bipolaronic) 
superconductor [1,6,17]. In case (iii), on the contrary, we find a situation which is similar 
to the BCS case: pairs of " c" electrons with opposite momenta and spins are exchanged via 
virtual transitions into local pair states. 

The generic phase diagrams for s-wave pairing symmetry plotted as a function of the 
position of the LP level A at fixed n are shown in Fig.l. In Fig. 2. we show the transition 
temperatures for the d x i _ y i -symmetry. 

In all the cases one observes a drop in the superfluid stiffness (and in the KT transition 
temperature) when the bosonic level reaches the bottom of the c-electron band and the 
system approaches the LP limit. In the opposite, BCS like limit, Tjf T approaches asymp- 
totically T^ FA , with a narrow fluctuation regime. Between the KT and MFA temperatures, 
phase fluctuation effects are important. In this regime a pseudogap in the c-electron spec- 
trum will develop and the normal state of LP and itinerant fermions can exhibit non-Fermi 
liquid properties [5]. 

A closer inspection of the Mixed-LP crossover indicates that when the LP level is lowered 
and reaches the bottom of the fermionic band an effective attraction between fermions 
becomes strong, since it varies as J 2 /(2A — 2/i) and ji pa A [12,13]. In this regime the 
density of c electrons is low and formation of bound c-electron pairs occurs. It gives rise 
to an energy gap in the single-electron spectrum independently of the pairing symmetry. 
We calculated the binding energies of c- electron pairs and found that T c Mi?A essentially 
scales with the half of their binding energy for A < 0. The superconducting transition 
temperature is here always much lower than the c-pair formation temperature (T C MFA ) and 
decreases rapidly with \A /D\. In such a case, the superconducting state can be formed by 
two types of coexisting bosons: preformed c- electron pairs and LP's [12,15]. 

Comparing T c vs A plots for various pairing symmetries one finds that in the case of nn 
hopping only, the d and s -wave pairings are favorable for higher concentration of c-electrons, 
while the s*-wave can be stable at low n c . The nnn hopping t 2 (with opposite sign to t) can 
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strongly enhance T c for <i-wave symmetry, moreover it favors the d and s-wave pairings for 
lower values of n c (compare Fig. 5) [12]. 

The region between T^ FA and T^ T ', where the system can exhibit a pseudogap, expands 
with increasing intersubsystem coupling \I \. As we have found [13], except for \Iq\/D <C 1 
the coupling dependences of T c Mi?A and Tf T are qualitatively different. T^ FA is an increas- 
ing function of |/ | for all the pairing symmetries. On the other hand, Tf T vs \I \ increases 
first, goes through a round maximum and then decreases (similarly as it is observed in the 
attractive Hubbard model). The position of the maximum corresponds to the intermediate 
values of \Iq\/D and it depends on the pairing symmetry as well as the values of A /D and 
n. For large |/ |, the Tf T are close to the upper bound for the phase ordering temperature 
which is given by 7rp s (0)/2. 

Concerning the evolution of superconducting properties with increasing n at fixed A 
one finds three possible types of density driven changeovers [12,8]: (i) for 2 > A /D > 0, 
"BCS" — ► Mixed — > "BCS"; (ii) for A /D > 2: "BCS" — ► "LP" and (iii) for A /D < 0: 
"LP" — ► "BCS". Only if the LP level is deeply located below the bottom of the c-band, 
the system remains in the LP regime for any n < 2. 

Let us also comment on the effects of a weak interlayer coupling on the calculated tran- 
sition temperatures [9,18]. In the KT theory the 2D correlation length behaves as follows 
for T > T KT ; £(T) = aexp (b/ ^Jt/T kt — lj , where b rs 1.5 and a is the size of the vortex 
core. If U c is the coupling energy per unit length between the planes and U c T KT , then 
the actual T c can be estimated by calculating the energy needed to destroy phase coherence 
between two regions of size ~ £ 2 in different planes i.e. T c ~ cU c (C,(T c )/a) 2 , where c is the 
interplanar distance. The resulting equations for T c can be solved asymptotically 

therefore T c is only weakly dependent on the interplanar distance c and is close to T KT , if 
U c Tkt- in the presence of the interplanar coupling there is no discontinuous jump in p s 
but a crossover from 2D like to 3D like (XY) behavior occurs. 
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IV. SUPERFLUID TRANSITION FROM THE PSEUDOGAP STATE 



Let us now consider the pseudogap behavior and present the recent evaluation of the 
superconducting transition temperature from a pseudogap state by going beyond the BCS- 
MFA. In our analysis we have applied a generalized T-matrix approach adapted to a two- 
component boson- fermion model [21]. Our approach is an extension of the pairing fluctuation 
theory of the BCS-Bose-Einstein crossover [19,20] developed previously for a one-component 
fermion systems with attractive interaction. The numerical results presented in Fig. 3. are for 
a 3D sc lattice assuming the tight-binding dispersion for fermions and bosons of the following 
form: e k = D(l -%), D = zt; J q = J 7q, Jo = zJ, % = [cos(k x ) + cos(k y ) + cos(k z )] /3, 
z = 6. 

The results are shown for both cases with and without the direct hopping of LP's J^. 
The calculated T c 's are much lower as compared to BCS-MFA results (these are given by 
Eq.(7)), and if J = 0, T c is strongly depressed as soon as the LP level is close to the bottom 
of the electronic band. In the pseudogap region the electronic spectrum is gapped, and the 
pseudogap parameter at T c for Ao > essentially measures a mean square amplitude of the 
pairing field (of the " c" electrons) . The values of pseudogap parameter at T c are comparable 
to the zero temperature gap values in the fermionic spectrum, except for the c-regime. 

With the direct LP hopping J /D = 0.1, which corresponds to m B = 10m p, the hard- 
core bosons can undergo a superfluid transition even without the intersubsystem coupling 
| Jo |- As we see from Fig. 3 in the presence of the boson- fermion coupling |/ | the transition 
temperature is enhanced in the mixed regime. 

In the self-consistent T-matrix approach the (amplitude) fluctuations of the order pa- 
rameter are included at the Gaussian level. Nevertheless, it is interesting to observe that the 
phase diagram for J = shown in Fig. 3 displays similar regimes as that of Fig.l determined 
in Sec. Ill from BCS and KT theories. 

V. FINAL REMARKS 

In conclusion, we summarize the important features of the model considered [11-13,5]. 
1. Well defined Bogoliubov quasiparticles can exist in the superconducting ground state. 
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However, above T c (in a mixed regime) local pairs coexist with itinerant fermions and 
the normal state properties deviate from Fermi liquid behavior. 

2. In the mixed regime, T C MFA < T < T^ T , the system will exhibit a pseudogap in the 
c-electron spectrum, which will evolve into a real gap as one moves to the LP regime. 
For A < 0, LP's coexist with preformed c-electron pairs, which have the binding 
energy E c h j2 oc Tf FA . 

3. For d-wave pairing, the superfluid density exhibits linear in T behavior at low T due 
to the presence of nodal quasiparticles. 

4. The Uemura-type plots i.e., the T c vs zero-temperature phase stiffness p s (0), are ob- 
tained for d, s* and s -wave symmetry in the KT scenario [12]. The reason for Uemura 
scaling T c ~ p s (0) is the separation of the energy scales for the pairing and for the 
phase coherence. [11,12]. 

5. The calculated T c 's in a 3D model beyond the BCS-MFA show crucial effects of pair 
fluctuations in the mixed and LP regimes. 

Some of our findings can be qualitatively related to experimental results for the cuprate 
HTS where a pseudogap exists. It has been suggested by ARPES experiments, that for 
underdoped cuprates the Fermi surface in the pseudogap phase is truncated around the 
corners due to the formation of preformed (bosonic) pairs with charge 2e, whereas the 
"electrons" on the diagonals remain unpaired [7a]. In the present two component model 
such a situation is obtained when LP's and c-electrons coexist in the mixed regime The 
linear T-dependence of the superfluid density has been observed experimentally in copper 
oxides and also in several organic superconductors. This points to an order parameter of 
rf a ,2_j / 2-wave symmetry and existence of nodal quasiparticles. In the present model the gap 
ratio is nonuniversal for all the pairing symmetries and can deviate strongly from BCS 
predictions (particularly in the d-wave case for which it is always enhanced) [12]. This 
feature is also found in several exotic superconductors. The Uemura plots and the scaling 
T c ~ Ps(0) reported for cuprates and organic superconductors can be reproduced within the 
model for extended s- and d-wave order parameter symmetry [11,12]. 
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FIG. 1. Phase diagrams of the induced pairing model as a function of Aq/D at fixed n derived for 
s-wave symmetry. Jq = and ti = 0. \Iq\/D = 0.25, D=4t. The dashed lines show the BCS-MFA 
transition temperature (upper for n = 1 and lower for n = 0.5), while the lines with circles and 
triangles show the KT transition temperatures calculated for n = 1 and n = 0.5, respectively. 
LPN-normal state of predominantly LP's, EM-electronic metal, LPS+ES-superconducting (SC) 
state, PG - pseudogap region. A weak interplanar coupling stabilizes the SC state. 
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FIG. 2. MFA and KT transition temperatures as a function of Aq/D at fixed n = 1 derived 
for d x 2_ y 2 - pairing symmetry. The dashed and dot-dashed lines show the BCS-MFA transition 
temperatures for t2 = and for t^lt = —0.45, respectively. The line with diamonds shows the corre- 
sponding KT transition temperature calculated for t2 = and the line with circles for ^/i = —0.45. 
\I \/D = 0.25, J = 0. 
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FIG. 3. Phase diagrams of the hard-core boson- fermion model as a function of Aq/D for s-wave 
pairing and sc lattice, n = 0.5, \Iq\/D = 0.5, D = 6t. The transition temperatures derived with 
a T-matrix approach are for two values of J, which are shown by the solid line (Jq/D = 0.1) 
and the line with symbols (Jo = 0), respectively. The dashed lines indicate BCS-MFA transition 
temperatures (upper for Jo/D = 0.1, lower for Jo = 0). LPN-normal state of predominantly LP's, 
EM-electronic metal, LPS+ES-superconducting (SC) state, PG - pseudogap region. 
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